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'We introduce a Green’s function method for handling radiative effects on false vacuum decay. 

In addition to the usual thin-wall approximation, we achieve further simplification by treating the 
bubble wall in the planar limit. As an application, we take the theory, extended with N 
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quadratic fluctuations about the classical soliton conhguration and the hrst correction to the soliton 
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I. INTRODUCTION 

The association of the still recently-discovered 125 GeV 
scalar particle [1, 2] with the Higgs boson of the Stan¬ 
dard Model (SM) places the stability of the electroweak 
vacuum under question [3-6]. This instability, arising at 
an energy scale around 10^^ GeV [7, 8], results from the 
renormalization-group (RG) running of the Higgs self¬ 
coupling, whose value is driven negative by contributions 
dominated by top-quark loops. State-of-the-art calcula¬ 
tions suggest that the electroweak vacuum is metastable, 
having a lifetime longer than the present age of the Uni¬ 
verse and lying at the edge of the stable region [7-10], 
where seemingly-small corrections may have a material 
impact upon predictions. The uncertainty in such pre¬ 
dictions remains, at present, dominated by that of the 
top-quark pole mass [11, 12]. Even so, it has been sug¬ 
gested [13-17] that, regardless of any improved precision 
in the experimental determination of the latter, the pres¬ 
ence of Planck-scale operators may weaken the claim of 
metastability. Nevertheless, having, as yet, no experi¬ 
mental evidence of additional stabilizing physics between 
the electroweak and Planck scales, it is provident to con¬ 
sider approaches to the calculation of tunneling rates that 
can consistently account for radiative corrections. 

The degree of vacuum metastability provides a strong 
criterion for the phenomenological viability of extensions 
to the SM. For example, supersymmetric scenarios can be 
ruled out if the electroweak symmetry-breaking vacuum 
decays into a color-breaking one in a timescale shorter 
than the age of the Universe [18-25]. In addition, tran¬ 
sitions between vacua can also occur at finite tempera¬ 
ture [26, 27]. In the context of early-Universe cosmol¬ 
ogy, this is of interest because the corresponding first- 
order phase transitions may leave behind relic gravita¬ 
tional waves [28-30]. Moreover, such phase transitions 
may turn out to be pivotal for generating the cosmic 
matter-antimatter asymmetry [31, 32]. As a consequence 
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of these applications and the wide range of phenomeno¬ 
logical models, there are now routine methods for com¬ 
puting transition rates at both vanishing and finite tem¬ 
perature [33, 34]. 

Vacuum transitions in scalar theories can be described 
in the following way [35-38]. In the event that there 
are two non-degenerate vacua, an initially homogeneous 
system lying in the false vacuum will spontaneously nu¬ 
cleate bubbles of true vacuum, leading to the produc¬ 
tion of domain walls or “kinks.” The latter are the 
topological solitons that interpolate between regions of 
true and false vacuum. The study of these “solitary 
wave” solutions to non-linear equations of motion (see 
e.g. Ref. [39]) has a long history [40-45], and archetypal 
examples of such field conhgurations arise in the sine- 
Gordon model [46, 47] and the theory with tachy- 
onic mass < 0. The semi-classical [38] and quan¬ 
tum [48] descriptions of false vacuum decay in the latter 
theory were presented in the seminal works by Coleman 
and Callan (see also Ref. [49]). Early expansion on these 
works included induced vacuum decay [50] and the incor¬ 
poration of gravity [51]. 

In order to decide whether a vacuum configuration is 
unstable, i.e. whether there exists a lowest-lying true vac¬ 
uum, it is often necessary to account for the impact of ra¬ 
diative corrections. This is of particular relevance when 
the appearance or disappearance of minima is entirely 
a radiative effect [52], such as occurs for the Coleman- 
Weinberg (CW) mechanism of spontaneous symmetry 
breaking [53] or in symmetry restoration at finite tem¬ 
perature [54-56]. In phenomenological studies, this is 
commonly done by calculating the tunneling rate from 
the effective potential [57, 58] of a homogeneous field 
configuration [34, 59, 60], which is subsequently pro¬ 
moted to a space-time-dependent configuration. This 
practice is problematic for two reasons. Firstly, the tem¬ 
poral and spatial inhomogeneity of the solitonic back¬ 
ground is not fully taken into account [61]. Secondly, in 
the presence of tachyonic instabilities, e.g. when there 
are non-convex regions in the tree-level potential, the 
perturbatively-calculated effective potential receives a 
seemingly-pathological imaginary part. The latter has 
been shown [62] to have a physical interpretation as a 
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decay rate for an initially homogeneous field configura¬ 
tion (see also Ref. [63]). However, this subtlety can be 
circumvented by using constructions such as the coarse¬ 
grained effective action [64-71] or, as is often employed in 
lattice simulations, the constraint effective potential [72] 
(cf. Ref. [73]). In addition, within the context of the SM, 
the standard RG improvement of the effective potential 
has recently been questioned [74]. For the above reasons, 
it is reasonable to conclude that the use of the effective 
potential to calculate transition rates is neither satisfac¬ 
tory nor justifiable and it is desirable to consider alterna¬ 
tive methods for determining quantum corrections, which 
can be applied to a wide range of models that feature 
vacuum decay. 

The first quantum corrections [48] to the tunneling 
rate are those arising from the functional determinant 
over the quadratic fluctuations about the classical soli- 
ton configuration. In the case of one-dimensional op¬ 
erators, these determinants may be calculated using the 
Gel’fand-Yaglom theorem [75], which may be generalized 
to higher dimensions in the case of radially-symmetric 
operators [76-78]. General numerical techniques may 
then be obtained [79, 80] for calculating tunneling rates 
beyond the so-called thin-wall approximation, in which 
the width of the bubble wall is much smaller than its 
radius. These approaches have also been applied to 
radially-separable Yang-Mills backgrounds [81, 82] and 
scenarios in curved spacetime [83]. Alternatively, as we 
will employ, the functional determinant may be calcu¬ 
lated by means of the so-called heat kernel method (see 
e.g. Refs. [84-87]), based upon the Schwinger proper-time 
representation and zeta function regularization [ 88 ]. Pre¬ 
viously, this approach has been used to derive approxi¬ 
mate analytic results for the one-loop fluctuation deter¬ 
minant beyond the thin-wall approximation [89], as well 
as one-loop corrections to sphaleron rates [90-92]. The 
latter have also been calculated by direct integration of 
the Green’s function [93-97]. 

Recently, it has been shown that properties of topo¬ 
logical solitons may be studied non-perturbatively using 
Monte Garlo and lattice simulations by considering cor¬ 
relation functions directly [98-100]. Other authors have 
proposed methods for calculating quantum corrections 
based upon functional renormalization techniques [ 101 ]. 

In this article, we derive an analytic result for the 
Green’s function of the theory in the background of 
the classical kink solution. Within the thin- and planar- 
wall approximations, we illustrate that this Green’s func¬ 
tion may be used to determine analytically the leading 
quantum corrections to both the semi-classical bounce 
action and the kink solution itself, thereby allowing us to 
compute the tunneling rate at the two-loop level, while 
isolating its diagrammatic interpretation. The latter cal¬ 
culation is performed within the context of a toy model 
extended with an additional N heavier scalars, where the 
parametric dependence on N allows the identification of 
a concrete example in which the calculated two-loop cor¬ 
rections dominate over the neglected higher-loop correc¬ 


tions. We illustrate that the problem of calculating these 
radiative corrections may be reduced to one of solving 
one-dimensional ordinary differential equations and inte¬ 
grals. Thus, we anticipate that this methodical develop¬ 
ment may have numerical applications in the study of the 
decay rates of radiatively-generated metastable vacua, 
such as occur in the massless GW model [53] or the Higgs 
potential of the SM. Similar methods based upon Green’s 
function techniques have been used previously to deter¬ 
mine self-consistent bounce solutions numerically in the 
Hartree approximation of the pure theory in both 
two and four dimensions [102-105]. 

The remainder of this article is organized as follows. 
In Sec. H, we review the calculation, a la Coleman 
and Callan [38, 48], of the classical “bounce” configu¬ 
ration, describing the semi-classical tunneling rate be¬ 
tween two quasi-degenerate vacua and its first quantum 
corrections. In Sec. HI, we outline a Green’s function 
method for the evaluation of the functional determi¬ 
nant over the quantum fluctuations about the classical 
bounce, making comparison with existing calculations. 
Subsequently, in Sec. IV, we illustrate that this Green’s 
function method may be used to calculate analytically 
and self-consistently the first quantum corrections to the 
bounce itself. In Sec. V, we conclude our discussions and 
highlight potential applications and future directions. Fi¬ 
nally, a number of mathematical appendices are included, 
outlining the technical details of the calculations summa¬ 
rized in Secs. HI and IV. 


II. SEMI-CLASSICAL BOUNCE 

We consider a real scalar field $ = ^(a;), with four¬ 
dimensional Euclidean Lagrangian £ = (9^$)^/2 -|- U 
and classical potential 

C/ = + + (1) 

The mass squared is m% = — < 0, 5 is of mass di¬ 

mension one, A is dimensionless, Uq is a constant and 
dfj_ = djdx^ denotes the derivative with respect to the 
Euclidean spacetime coordinate = (x, 0 : 4 ). Throngh- 
out, we omit spacetime and field arguments for notational 
convenience when no ambiguity results. 

The classical potential in Eq. (1) has non-degenerate 
minima at 

ip = v± = -w, (2) 

as depicted in Fig. 1 (left panel), where we have defined 
V = y/OpPjX and v = {‘ig)/{2X). The separation of the 
minima Au = — V- = 2d and the difference in their 

energy densities AU = — C/„_ = 2s may be written 
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U{ip) -U{^) 




comes an 5'0(1, 3) symmetry, with the bubble expanding 
along the hyperbolic trajectory — c^t^. 

The bounce action is 


B 



1 

2 



+ U{ip) \ , 


( 6 ) 


which can be written as B = -Bsurface + ^vacuum, where 


FIG. 1. The classical potential U (left panel) and the inverted 
potential — U (right panel). The arrow (right panel) indicates 
the trajectory of the “bounce” in imaginary time r. 


in terms of the parameters 


d = V 



V , 


e 

d 



gv^ 
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(3a) 

(3b) 


where the approximations are valid in the limit v ^ v, 
i.e. 8A/3. For g —>• 0, e: —)■ 0 and the minima at 

If = ±v become degenerate, as we would expect. 

Finally, the constant Uq is chosen so that the potential 
vanishes in the false vacuum at ip = +v, requiring Uq = 
— gv^ jd) and thus giving the barrier height to be 
h = Uq + 2£ K, (^u/2)^ + e and U{—v) = — gv^/3. 

The semi-classical probability for tunneling between 
the false (tp = -I- u) and true {tp = —v) vacua and its first 
quantum corrections were described in the seminal works 
by Coleman and Callan [38, 48], which we now review. 
The classical equation of motion 


V 


.^surface — ^TT R 


2 c>3 


dip 
dip — 
dr 


Sv 


fR 

= 27r^ / drr^U{—v) 

Jo 


(7a) 

(7b) 


are the contributions from the surface tension of the bub¬ 
ble and the energy of the true vacuum, respectively. In 
writing Eq. (7a), we have used the fact that for df^ip ^ 0, 
i.e. for r ^ R, we may show that the bounce ip satisfies 
the virial theorem 
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i.e. it is the configuration of zero total energy density. 
Notice that there is no contribution to the bounce action 
[Eq. (6)] from the exterior of the bubble, since the choice 
of the potential Eq. (1), viz. Uq, ensures that the false 
vacuum has zero energy density. 

In the thin-wall approximation, we may safely neglect 
the damping term in Eq. (5) and the contribution from 
the cubic self-interaction gip^, as will be the case for the 
remainder of this article. We then obtain the well-known 
kink solution [40] 


-d‘^ip + U'{ip) = 0 , (4) 


ip{r) = vtanh ^{r — R)] , (9) 


where ' denotes the derivative with respect to the field 
ip, is analogous to that of a particle moving in a po¬ 
tential — U. The boundary conditions of the “bounce” 
are ip\xi^±oo = +u and ip\xi=o = 0, where ' denotes 
the derivative with respect to Xj^. These correspond to 
a particle initially at -I- u rolling through the valley in 
— U, reaching a turning point close to — v, before rolling 
back to -bv, see Fig. 1 (right panel). Finally, in order to 
ensure that the action of the bounce is finite, we require 
T’IixKoo = +^'■ 

Translating to four-dimensional hyperspherical coordi¬ 
nates, Eq. (4) takes the form 


d 2 

dr 2 


3 d 
r dr 


ip + U'{ip) 


0 , 


( 5 ) 


where + xf. The boundary conditions be¬ 

come ip\r^oo = + v and dip/dr\x=Q = 0, where the lat¬ 
ter ensures that the solution is well defined at the ori¬ 
gin. Thus, the bounce corresponds to a four-dimensional 
bubble of some radius R, which separates the false vac¬ 
uum {ip = -b u) outside from the true vacuum inside 
{ip = —v). Analytically continuing to Minkowski space- 
time (x 4 = ixo), the 0(4) symmetry of the bounce be- 


with 7 = The radius R of the bubble is then ob¬ 

tained by extremizing the bounce action Eq. (6), that is 
by minimizing the energy difference between the surface 
tension of the bubble and the true vacuum. This gives 


R = 


127 

gv 


( 10 ) 


By considering the invariance of the bounce action 
in Eq. (6) under general coordinate transformations, 
i.e. ip ^ ip + x^df^ip, we may show that 


B = -tt^R^ 
2 


' — V 


dr 


This is to say that 
B. 


vacuum 


= -^S.urface, 


( 11 ) 


( 12 ) 


in which case we find 

I TT^ 

B = - -^vacuum = - — R'^ U{-v) 
3 6 


A ■ 

(13) 



4 


The decay rate of the false vacuum, i.e. the probabil¬ 
ity per unit time for the nucleation of a bubble of true 
vacuum, has the generic form [38, 48] 

r = AVe-^/^. (14) 

Here, V is the three-volume within which the bounce 
may occur, arising from integrating over the center of 
the bounce, and A contains the quantum corrections to 
the classical bounce action B that are the subject of the 
remainder of this article. 

The tunneling probability in Eq. (14) may be obtained 
from the path integral 

Z[0] = y[d$] , (15) 


via 


r = 2|ImZ[0]|/T, (16) 


where T is the Euclidean time of the bounce. 

In order to evaluate the functional integral over $, we 
first expand around the classical bounce (/?, whose equa¬ 
tion of motion (4) is obtained from 


( 54 ) 


^ 


0 . 


(17) 


Writing $ = where the factor of is written 

explicitly for bookkeeping purposes, we find 


S'[4)] = S'[<^] -I- ^y'd'^a; <i)(x) G ^ {ip-.x) ^{x) A , 

(18) 

where S'[v3] = B \s the classical bounce action and 


G '^{ip]x) 


54>2(a;) 




-AW+[/"((/); a;) , (19) 


in which is the four-dimensional Laplacian. 

Before proceeding to perform the functional integra¬ 
tion over the quadratic fluctuations about the bounce, 
we must consider the spectrum of the operator G“^(iy9; x), 
which is not positive dehnite. By differentiating the equa¬ 
tion of motion Eq. (4) with respect to x^ and comparing 
with the eigenvalue equation 

( _ AW + C/"(<p)) = \{n}hn} ( (20) 


it is straightforward to show that there exist four zero 
eigenmodes 4>^ = transforming as a vector of 

S'0(4) and resulting from the translational invariance of 
the bounce. The normalization M follows from Eq. (11), 
since 



( 21 ) 


Differentiating Eq. (5) with respect to r and subse¬ 
quently setting r = i? in those terms originating from 
the damping term, we can show that there also exists 
a discrete eigenmode (/)o = B~^^‘^dr(p. This eigenmode 
transforms as a scalar of S'0(4), corresponding to dilata¬ 
tions of the classical bounce solution, and has the nega¬ 
tive eigenvalue 


A = = - A f22l 

BSR^ ■ ( ) 

It is this lowest mode that is responsible for the path in¬ 
tegral in Eq. (15) obtaining the non-zero imaginary part 
in Eq. (16) [106]. 

Alternatively, we may solve the eigenvalue problem di¬ 
rectly in hyperspherical coordinates (see Appendix B), 
by making the substitution (^{„} = Neglecting 

the damping term and setting r = R in the centrifugal 
potential, we obtain the eigenspectrum 

— 7 (4 - n ) -I-- . (23) 

The radial parts of the eigenfunctions are the associ¬ 

ated Legendre polynomials of the Hrst kind and of order 
2, i.e. Pi^^ip/v). Thus, demanding normalizability, the 
quantum number n is restricted to the set { 1 , 2 }. 

From Eq. (23), we see that the negative mode corre¬ 
sponds to Ao = A 20 (n = 2, j = 0) and the zero modes 
correspond to A 21 (n = 2 , j = 1 ), having degeneracy 
(j + 1)^ = 4. The lowest two positive-definite eigenvalues 
are Aio = 272 — 3/i?^ (n = 1 , j = 0 ) and An = 2j^ 

(n = 1, j = 1). Thus, for R large, the “continuum” 

of positive-definite modes begins at Aio « An = 272 , 
cf. Ref. [ 86 ]. 

In order to perform the functional integral over the 
five ne|ative-semi-deHnite discrete modes, we expand 
$ = where 4)+ comprises the continuum 

of positive-definite eigenmodes. The functional measure 
then becomes 

4 

[d4>] = [d(/)+] J](2^fi)-i/2da, . (24) 

i=0 

The functional integral over the four zero eigenmodes 
{i = 1,... ,A) is traded for an integral over the collective 
coordinates of the bounce [107] (see Appendix A) and 
yields a factor 


The integral over the negative eigenmode (i = 0) may 
be performed using the method of steepest descent, giv¬ 
ing an overall factor of —t|Ao|“^^^/2. Here, the overall 


^ We note that this substitution differs from that used in Ref. [48]. 


Thus, (t)y. = B 
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sign is unphysical [48] and depends on the choice of an¬ 
alytic continuation, thereby justifying the modulus sign 
in Eq. (16). 

Finally, the Gaussian integral over the continuum of 
positive eigenmodes may be performed in the usual 
manner and we obtain 

iZ\0] = Apde#^ 

(26) 

in which det*'^^ cf. Ref. [86], denotes the determinant 
calculated only over the continuum of positive-definite 
eigenmodes, i.e. omitting the zero and negative eigen¬ 
modes, whose contributions are included explicitly. In 
addition, we have normalized the determinant to that 
of the operator G~^{v), evaluated in the false vacuum. 
Substituting Eq. (26) into Eq. (16), we find the tunneling 
rate per unit volume 

(27) 

where 

= itr(5)(lnG-i((p)-lnG-i(u)) , (28) 

contains the one-loop corrections from the quadratic fluc¬ 
tuations around the classical bounce. Here, tr^^^ indi¬ 
cates that we are to trace over only the positive-definite 
eigenmodes. 

III. GREEN’S FUNCTION METHOD 

In this section, we outline the derivation of the Green’s 
function of the operator in Eq. (19). The technical de¬ 
tails are included for completeness in Appendix B. Sub¬ 
sequently, we use this Green’s function to evaluate the 
functional determinant in Eq. (26) and obtain the correc¬ 
tion from quadratic fluctuations. In addition, we calcu¬ 
late analytically the tadpole contribution to the effective 
equation of motion and point out that this may be used 
to calculate the first quantum corrections to the bounce. 
We have the inhomogeneous Klein-Gordon equation 

+U''{(p;x))G(ip;x,x') = S^'^^x — x') , (29) 

where — x') is the four-dimensional Dirac delta 

function. Working in hyperspherical coordinates and 
writing xj^''^ = where e,,(,) are four-dimensional 

unit vectors, the Green’s function may be expanded as 

- OO 

G{ip;x,x') = ^'^{j + l)Gj{Lp-,r,r')Uj{cos0) , (30) 

where cos6l = • e^' and Uj{z) are the Chebyshev poly¬ 

nomials of the second kind (see Appendix B). The radial 


functions Gj{r,r') satisfy the inhomogeneous equation 


_ d_d_ j(j + 2) 

dj.2 ^ j,2 


+ [/"(r) G,{^-ry) = 


6{r — r') 


For the thin wall, we safely neglect the damping term 
and approximate the centrifugal term by j(j-|-2)/i?^. For 
self-consistency of this approximation, we also replace the 
discontinuity on the rhs of Eq. (31) with 5{r—r')/E?. For 
generality of notation in what follows, it is then conve¬ 
nient to define 

G(u,u',m) = R^Gj{ip;r,r') , (32) 

being a function only of the normalized bounce 

m(') = ^ = tanh [ 7 (r^'^ — R)] (33) 

V 

and the parameter 


m = 2(1-1- 


J'(J + 2 )V 


V 472772 j 

The full Green’s function may then be written 
G{ip]x,x') = G(u,u',9) 


^ 00 

+ ^)Uj{cos9)G{u,u',m) 


27r2i?3 


With the above approximations, the Ihs of Eq. (31) is 
of Poschl-Teller form [108], having general solutions that 
may be expressed in terms of the associated Legendre 
functions (see Appendix B). We are then able to find the 
full analytic solution 


G{u,u',m) = ^ d{u — u') 

2'ym [ 




X 1-3 


-I- (m o u') , (36) 


x(i-3 <(7">3+’" + “) ) 

\ (1-I-m)(2-I-m) ) 

( (1 — m')(1 — m-I- m') \ , , , 

X ( 1 — 3 -—-r—-r- j + {u ^ u) , (36) 

where d(z) is the generalized unit-step function. 

Taking the coincidence limit u = u', 9 = 0, the local 
contribution to the Green’s function G{u) = G(u, u, 0) in 
Eq. (35) takes the form 

^ 00 

G{u) = 2 ^ 2 + ^fG{u,m) , (37) 


where 

G{u,m) = G(u,u,m) 


1 + 3(1_„2)^ 


(—l)”(n — 1 — v?) 
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FIG. 2. Schematic representation of the planar-wall approx¬ 
imation, illustrating the alignment of the three-dimensional 
hypersurface and associated coordinate system. 


In Eq. (38), the summation over n = 1,2 corresponds to 
the contributions from the two towers of positive-definite 
eigenmodes of the operator G~^{(p-,x), see Eq. (23). 

For R large, we may approximate the summation over 
j by an integral over a continuous variable k ~ (see 
Appendix B). In which case, we obtain 

1 

dkk^G{u,m), (39) 

27r Jq 


with 


A. Quantum-corrected bounce 

Before making use of the Green’s function calculated in 
the preceding section, we first derive the equation of mo¬ 
tion for the quantum-corrected bounce. This calculation 
was first suggested by Goldstone and Jackiw [42] and, in 
the following sections, we will illustrate that, within the 
thin- and planar-wall approximations, it may be com¬ 
pleted analytically. 

The one-particle irreducible (IPI) effective action [57] 
is given by the Legendre transform 

r[0] = — J di^x J{x)(j){x) , (44) 


where 


^(x) 


^ 6 In Z[J] 
^ SJix) 


is a functional of the source J{x) = 


and 


(45) 


Z[J] = J[d^] exp — — Jd'^x J{x)^{x) 


(46) 

In order to obtain the quantum corrections to the 
bounce tp, we wish to evaluate the functional integral 
in Eq. (46) by expanding around the configuration 
which is the solution to the quantum equation of motion 


’" = 2([l + Ay. (40) 

The continuum limit described above is entirely 
equivalent to the so-called planar-wall approximation. 
Therein, for R large, we align a set of coordinates (^_l, Zj[) 
with the bubble wall, as shown in Fig. 2. We may then 
Fourier-transform with respect to the coordinates zy that 
lie within the three-dimensional wall, introducing a three- 
momentum k, i.e. 

/ d^k / 

—— G{(p;z,z',k) , (41) 


sm 

5(j) 




0 . 


(47) 


Here, the superscript “(1)” indicates that contains 
the first quantum corrections to p. It follows from 
Eq. (47) that cannot extremize the classical action 
in the absence of the source J, i.e. 


dS[^] 

(5$(a;) 




J{x) ^ 0. 


(48) 


Writing $ = where the factor of 

is again written explicitly for bookkeeping, we proceed 
as in Sec. II, expanding 


where we have let z = z± for notational conve¬ 
nience. The three-momentum-dependent Green’s func¬ 
tion G{p;z,z',k) satisfies the inhomogeneous Klein- 
Gordon equation 


S'[4>] = S[p^^^]+h^/^Jd^x J{x)¥^'>{x) 

+ ^ Jd‘^x^^^\x)G-\p''^'>;x)^^^\x) + -- - , (49) 


(^ —8“^ + k'^ + U"{p; z))G{p-, z, z'jk) = S{z — z') . (42) where 


We may then show straightforwardly that 

G{p;z,z',k) = G{u,u,m), (43) 

where G{u,u',m) is as defined in Eq. (36) with m given 
by Eq. (40). This planar-wall approximation is employed 
for the remainder of this article. 


G-\p^^^;x) 


52A[$] 

S^^lx) 


«J) = (^( 1 ) 




(50) 

We now write = p + h6p and expand about 

the quadratic fluctuations evaluated around the classi¬ 
cal bounce p. Thus, in performing the functional inte¬ 
gral, we consider the same spectrum of negative and zero 








7 


eigenmodes as in Sec. II. Finally, by expanding the effec¬ 
tive action in Eq. (44) around = 4> ~ (see 
Ref. [109]), we obtain 


r[(^W] = ^[^( 1 )] + *^ + 


+ - 1 „ 


Ao det(^) 


{VTfi^) G-^v) 


where 


i y d^x S(p{x) 


In 


det^^^G 


S^W{x) det(^)G-i(^) 


(51) 


(52) 


5 ( 1 ) =(/ 


Functionally differentiating Eq. (51) with respect to 
we obtain the equation of motion for the corrected 
bounce 


Hence, we obtain the tunneling rate per unit volume 

r/V = 2|Ime"^[^‘'’'/'^|/(yT) 


= f-^) 

\ 2 Trh I 


X exp 


(27)'|Ao|-^ 


1 


- - {3 + 
h 


(61) 


where B is the classical bounce action; given in 

Eq. (28), contains the corrections from quadratic fluc- 
tations about the classical bounce; and B^‘^\ given in 
Eq. (60), contains the contribution arising from the quan¬ 
tum corrections to the bounce itself. We note that 

5 ( 2 )' = -2R(2) ^ (02) 

such that the 0{h) corrections to the quadratic fluctua¬ 
tions flip the sign of the contribution to the bounce action 
from the 0{h) corrections to the bounce itself. 


-d‘^(p^^\x)+ U^g((p^^^;x) = 0, (53) 

where 

= U'((p^^^;x)+hn{(p;x)(p{x) , (54) 

containing the tadpole contribution 


B. Tadpole contribution 

We will now proceed to calculate explicitly the tadpole 
contribution appearing in Eq. (55). 

Introducing an ultraviolet cutoff A, the k integral can 
be performed in Eq. (39), and we obtain 


n(7’;a;) = ^G((p;x,x) . (55) 

Comparing the functional derivative of Eq. (51) with 
Eqs. (47) and (48), we see that this evaluation of the 
effective action is self-consistent so long as the source 

J(x) = — hll{ip;x)ip{x) , (56) 

which is, as expected, non-vanishing. 

We may show that the correction to the classical 
bounce Sip satisfies the equation of motion 

G~^{ip;x) 6ip{x) = —B{ip-,x)ip{x) . (57) 

The corrected bounce action contains contribu¬ 

tions at order . Specifically, 

^[<^( 1 )] = 

f 

-I- — d‘^xS(f{x)G~^{(p;x)Sip{x) + 0{h^), (58) 

where we have used Eqs. (17) and (19). Thus, using 
Eq. (57), we may write 

S'[(^(i)] = B + h^B^'^\ (59) 

where 

(p{x)Il{if, x)5ip{x) . (60) 


Giu) = 


v2 r \2 ^,2 

^+2-(l-3u^)ln 


87r2 


T 

A 2 


— •\/37rM^(l — 


(63) 


We choose to define the physical mass and coupling in 
the homogeneous non-solitonic background.^ The renor¬ 
malization conditions are then as follows: 


d^Ueei^p) 

dip^ 


2,^2 

= + 2 - = 


d^Ueffip) 


Lp—V 


A, 


, 


(64a) 

(64b) 


where UeS is the CW effective potential [53]. The result¬ 
ing mass and coupling counterterms are 


6m^ = — 


SX = - 


167r2 

3A2 

327r2 



(65) 

>“^+ 5 ) ■ 

(66) 


2 It is natural to define the renormalized quantities in the false 
vacuum, since this is where the physical measurements of these 
quantities are performed. If it were the case that such measure¬ 
ments were taking place in the true vacuum, or indeed within 
the wall itself, then the decay rate would be of little concern. 



We then arrive at the renormalized tadpole correction 


n^(u) = ^G{u) + dw?' 


2 

3A7^ 

IGtt^ 


6 + (1-u2)(5 - 


(67) 


C. Functional determinant 


is the inverse Laplace transform with respect to s, with 
C indicating the Bromwich contour. 

We may perform the integrals in Eq. (74) analytically, 
proceeding in order from right to left and beginning with 
the inverse Laplace transform. We then obtain the un¬ 
renormalized correction to the bounce action 


^(1) 




A2 



(77) 


We may calculate the traces appearing in the exponent 
of Eq. (27), which arise from the functional determinant 
of the operator in Eq. (26), by using the heat 

kernel method (see e.g. Ref. [87]). Specifically, the trace 
may be written in the form 

/ poo 1 

d^x / — K{(P]x,x\t) . 

( 68 ) 

The heat kernel K{Lp-, x, x'\ r) is the solution to the heat- 
flow equation 

X, x'l r) = G~^{ip]x)K{(p]x^x'\T) (69) 

and satisfies the condition K{ip]x,x'\ 0) = S^^\x — x'). 

It is convenient to work in terms of the Laplace trans¬ 
form of the heat kernel 

poo 

JC{(f]x,x'\s) = / dr Rr(V3; X, x'l r) , (70) 

Jo 

which is the solution to 

( — d"^ + s + U"{(p;x))lC{(p;x,x'\s) = 6^'^\x — x'). (71) 
In the planar-wall approximation, we take 

IC{(p;x,x'\s) = e*‘"'^^ii"^ii^/C((p;z,z',k|s) , 

(72) 

where IC{(p; z, z\ k| s) satisfies 

{-dl + k"^ + s + U''{(p\z))lC{ip]z,z','k\s) = 5{z - z') . 

(73) 

Comparing Eq. (73) with Eq. (42), we see that 
K,{(p] z, z', k| s) is nothing other than the Green’s function 
G{u, u', m) in Eq. (36) with the replacement —)■ -|- s 
in m, see Eq. (40). Thus, we may write 



where we have defined 


G{u,m) = G{u,m) — G{l,m) (75) 

and 

C~^[Giu,m)]{T) = e"®^G(M,m) (76) 


The technical details of the relevant integrations are in¬ 
cluded in Appendix B. Adding the counterterm 

= J d'^x — v^) + ^6X(^ip‘^ — 


= B 


we obtain the final renormalized result 


= -B 




VibTT^y V3v^ 


21 


(78) 


(79) 


In Appendix B, we reproduce this result by the method 
presented in Ref. [86].^ 


IV. RADIATIVE CORRECTIONS TO THE 
BOUNCE 

We now discuss an example of the role played by loop 
corrections to the bounce itself. Within the perturba¬ 
tion expansion, one should expect that these lead to 
second-order corrections to the classical action of the 
soliton simply because the latter is evaluated for a sta¬ 
tionary path. There are, however, important situations, 
in which all one-loop contributions must be resummed 
in order to capture the leading quantum corrections 
to the action. Examples include situations where the 
symmetry-breaking minima of the potential emerge ra- 
diatively through the CW mechanism [53]. In the ab¬ 
sence of a soliton, this implies that the classical solution, 
i.e. the homogeneous expectation value of the field, has to 
be found consistently by minimizing the one-loop effec¬ 
tive potential as a function of the field expectation value 
itself. Analogously, in order to find the decay rate of the 
false vacuum, the bounce must be computed consistently 
from the one-loop effective action, which is a functional of 
the bounce itself. The methods presented in this article 


® Using the same renormalization conditions as in Eq. (64), 
Ref. [86] finds (in the notation employed here) 


= -B 





Repeating the analysis presented therein, as outlined in Ap¬ 
pendix B, we find a result in agreement with Eq. (79) reported 
here, suggesting a numerical error in the factor of 50 above. 
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reduce the problem of tunneling in radiatively-generated 
potentials to one-dimensional ordinary differential equa¬ 
tions and integrals. It is anticipated that it should be 
possible to derive numerical solutions in future work. 

For the purpose of illustration, however, we remain 
herein on the ground of analytic and perturbative ap¬ 
proximations. In order to enhance the corrections to the 
bounce compared to other quantum effects that appear 
at second order in perturbation theory, we extend the 
model in Eq. (I) with N copies of an additional scalar 
field X by adding to the Lagrangian the terms 

I • ( 80 ) 

Here, we have chosen the coupling A to be identical to the 
self-coupling of $ for the sake of simplicity in the Green’s 
function of the x fields. Since (xi) = 0, the additional 
scalars do not impact upon the classical bounce in Sec. II 
or the discussion of the Green’s function in Sec. III. 

The Klein-Gordon equation for Xi takes the form 

q2 I 2 I 2 

-d^ +m^ + 

Gomparing with that of $ in Eq. (29), we see that the 
Green’s function S{u,u',m) may be obtained straight¬ 
forwardly from G{u,u',m) in Eq. (36) by making the 
replacement 


S{(p;x,x') = 6^‘^\x — x'). (81) 


m —>■ VG 



6-^2 


1 

2 


(82) 


The renormalized tadpole contribution from each Xi 
field, integrated over the three-momentum k, is given by 


where we have assumed ^ 7 ^ for simplicity. The full 
form of S{u) and the relevant counterterms are provided 
in Appendix G. 

The renormalized correction to the classical bounce 5(p 
is governed by the equation of motion 


dr^ 



(n^(n) + AfE^(u))¥>, (84) 


cf. Eq. (57). We obtain the solution by making use of the 
Green’s function G(u,u',2) = G(u,M',TO)|fe = 0 : writing 


6(p{u) 


V , u'G{u, u', 2) 

iJJ“ 




where we have used Eq. (33) in order to substitute ip. 

We note at this point that G{u,u',m) is singular as 
k ^ 0 (or, equivalently, m —)■ 2). Nonetheless, the 



K(r-R) 

FIG. 3. The correction to the bounce Sp as a function of 
7 (r — R) for jm^ = 0 (solid), 0.5 (dashed), 1 (dash- 
dotted) and 1.5 (dotted). 



K{r-R) 



FIG. 4. The corrected bounce p>-\-5ip as a. function of 7 (r — i?) 
for lm\ = 0 (solid), 0.5 (dashed), 1 (dash-dotted) and 
1.5 (dotted). We see clearly that the impact of the tadpole 
correction is to broaden the bubble wall. 


integral in Eq. (85) remains finite, since G{u,u',m) is 
multiplied with an odd function, whereas the singularity 
resides in its even part. It is therefore useful to define 

G°'^‘^(u,u') = ^(G{u,u',2)-Giu,-u',2)y (86) 

Within the domain 0 < u, u' < 1, this function can be 
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expressed as 


= 'd{u-u') 


1 1-u^ 




1 — It' 


327 1 - w'2 

+ (it O u) . 


2u'(5-3u'2) 

(87) 


Defining in addition 


Poiu) =7 / du' ^ 
7-1 1 - 


'-1 

1 - tt^ 


2m 1 + m 

-^ +ln- 

1 — m2 1 — m 


Pi(m) = dM' m'G°‘^‘^(m,m') 


1 — M^ , 1 + M 


In 


1 — M ’ 


P2{u) ~ "fj dM' m'^G°‘^‘^(m, m') 


1-1 
1 - M^ 


, 1 + M 4 

In- u 

1-M 3 


we find the result 
3Av 


5ip{u) = — 
+ 5 


167r2 

167 ^ 


^7V + 1 )po(m) 

m ^ ' 


3m2 


( 88 a) 


( 88 b) 


(88c) 


iV + 1 )pi(m) - ( ^ + -^']p2iu) 


mt V3^ 


(89) 


In Fig. 3, we plot 6ip as a function of 7 (r — R) for a range 
of values of We see from Fig. 4, which plots 

the corrected bounce (p + Sip for the same range, that 
the impact of this correction is to lower the height and 
broaden the width of the bubble wall. We note that this 
behaviour is in qualitative agreement with the results of 
the self-consistent numerical analysis in Ref. [102], pre¬ 
sented there for the pure theory in 1 -|- 1 dimensions. 

Substituting Eq. (89) into Eq. (60), we find the correc¬ 
tion to the bounce action 


^(2) 


''d^x ip{u) (n^(M) -b 1 VE^(m)) 6ip{u) 


B 

J 



291 37 TT 

2897 TT \ 72 
42 ^3 / m2 


5 TT^ 

M y 

N 


5829 7 ^ 
14 



(90) 


In order to obtain a finite result for Eq. (90), we have 
added to Go the correction 


6U0 


9 

4 



2 2 
7 V 



N + 1 


2 


(91) 




FIG. 5. Diagrammatic representation of various contribu¬ 
tions to the effective action: (A) is the one-loop term 
(B) is the O(X^N^) contribution to and (C) and (D) 

are 0{S?N) terms. Solid lines represent the propagator 
G{ip\x,x'), dotted lines S{p\x,x'). Crosses denote insertions 
of the bounce p. 


ensuring that the potential continues to vanish in the 
false vacuum. 

The corrections appearing in Eq. (90) should be com¬ 
pared to the renormalized logarithm of the determinants 
of the Klein-Gordon operators of the y fields in the back¬ 
ground given by (/?, which are given by 


^167r2y 15 






N. (92) 


In comparison, the leading term in Eq. (90) is suppressed 
by a factor ~ A^^/m^/( 167 r^). The one-loop corrections 

and are both negative, thereby increasing the 
tunneling rate. It is interesting to note that, although 
the contribution to the tunneling action from the 
corrections to the bounce itself is positive, the net con¬ 
tribution of B^^'> is still negative, again increasing 

the tunneling rate. 

In Fig. 5, we present a diagrammatic representation 
of the corrections to the bounce action. It is also use¬ 
ful in order to see that there appear no contributions 
of 0{X^N^) relative to the bounce action B in addition 
to those from In order to avoid proliferation, we 

only show the leading contributions in 1 /A^ for a given 
type of diagram. At one-loop order, there is the vac¬ 
uum bubble in terms of the propagator S of the y fields. 
Fig. 5(A), which gives the contribution 0{XN) relative 
to B from B^'^ in Eq. (92). On substituting Sip in the 
form of Eq. (85) into the action [Eq. (58)], we see that the 
diagram corresponding to the 0{X^N^) term in B^^^/B 
is given by Fig. 5(B), where, when counting the powers 
of A, one should note that each explicit insertion of ip 
contributes a factor of \/VX. Finally, at two-loop or¬ 
der, there are the diagrams Figs. 5(C) and (D), which we 
do not compute but yield contributions of 0{X^N) rela- 
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live to B. These contributions are therefore suppressed 
by a relative factor of 1/A^ relative to the 0 {}?N‘^) in 
as is familiar from the standard approximation 
scheme known as the 1/iV expansion [110]. We should 
remark that these arguments do not hold, of course, for 
the contribution to from the $ tadpole, which we 
include here for completeness. The latter is formally the 
same order as other two-loop diagrams, involving only 
$, that are not captured in the IPI approximation em¬ 
ployed here. This observation is true also of the Hartree 
approximation for the pure theory analyzed numer¬ 
ically in Refs. [102-105]. Nevertheless, these additional 
two-loop diagrams remain subdominant compared to the 
0 {\^N) and 0 {X‘^N‘^) contributions from the y tadpole 
in Eq. (90). 

Finally, we note that approximating Sip as a small 
perturbation to tp , using Eq. (85), requires for consis¬ 
tency that 6A^A7^/(m^7r^) <C 1, such that within the 
range of validity of present approximations, we cannot 
obtain Nevertheless for large N, 

S(2) +S(2)' 

can be the dominant two-loop contribution 
to the effective action. 


V. CONCLUSIONS 

Within the context of A$^ theory, we have described 
a Green’s function method for handling radiative effects 
on false vacuum decay. By this means and employing 
the thin- and planar-wall approximations, we have been 
able to calculate analytically and in a straightforward 
manner both the functional determinant of the quadratic 
fluctuations about the classical soliton configuration and 
the first correction to the conhguration itself. 

This Green’s function method is well suited to numer¬ 
ical evaluation and, as a consequence, should be appli¬ 
cable to potentials of more general form. As such, we 
anticipate that it may be of particular use when the 
non-degeneracy of minima is purely radiatively gener¬ 
ated. Examples of the latter include the spontaneous 
symmetry breaking of the massless CW model [53] or 
the instability of the electroweak vacuum. Other appli¬ 
cations might include the calculation of corrections to in¬ 
flationary potentials in the time-dependent inflaton back¬ 
ground, for instance in inflection-point or A-term infla¬ 
tion [111-114], which exploit the flat directions and sad¬ 
dle points of the MSSM potential. Furthermore, the use 
of Green’s functions naturally admits the introduction 
of finite-temperature effects or extension to non-trivial 
background spacetimes. 

Green’s functions have proved to be central objects 
within perturbative calculations throughout quantum 
field theory and it is therefore unsurprising that we find 
these suitable to treat solitons in A$^ theory as well. We 
take this as an encouragement that further theoretically- 
and phenomenologically-interesting systematic results on 
false vacuum decay may be within reach. 
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Appendix A: Zero-mode functional measure 

In order to perform the functional integration over the 
zero modes, we insert four copies of unity in Faddeev- 
Popov form [107]: 

1 = • (Al) 

Here, y is not summed over and 
/(i/m) = Jd^x ^{x-y)dl^^ip{x-y) = , (A2) 

where we recall that 

4 

^ + + . (A3) 

It follows that 

= -Jd'^x{dl^'>ipix-y)Y = -B, (A4) 
ignoring terms that are formally Thus, 

\ = B^Ay^ S{B^l^ay) = B^/^ Jdy^ S{a^) . (A5) 

We then have 



Appendix B: Green’s function 

In this appendix, we include the technical details of the 
calculations outlined in Secs. Ill and IV. All functional 
identities used in what follows may be found in Ref. [115]. 
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1. Expansion in hyperspherical harmonics 

In d dimensions, the Green’s function satisfies the in¬ 
homogeneous Klein-Gordon equation 

= 5^'^\x-x'), (Bl) 

where {x — x') is the Dirac delta function and is 
the Laplacian. Given the 0{d) invariance of the bounce 
ip, it is convenient to work in hyperspherical coordinates, 
in which case the Laplacian takes the form 

A(d) = + Agd-i , (B2) 

where Agd-i is the Laplace-Beltrami operator on the d—1 
sphere. 

We proceed by performing a partial-wave decomposi¬ 
tion of the Green’s function: 

G^‘^'>{ip-,x,x') = ^Gj(v3;r,r')K/{^}(er')^jW(er) , 

(B3) 

where x = re^, x' = r'e^' and Yj^iy{er) are the hyper¬ 
spherical harmonics (see e.g. Ref. [116]), satisfying the 
eigenvalue equation 

Asd-iYj{iy = — j{j + d — 2)Yj^iy , (B4) 

with {G\ = i\, £ 2 , ■■■, £d- 2 - The hyperradial function 
Gj{(p;r,r') satisfies 


For d = 2, K = 0, and we have 

. (BIO) 

2 

where Tj (z) is the Chebyshev polynomial of the first kind. 
We then obtain 

- 00 

G’''^\(f,x,x') = -'^cosj9Gj{(f,r,r') , (Bll) 

TT ^' 

1=0 

where we have used the trigonometric form Tj{cos6) = 
cos jd. 

For d = 3, K = 1/2, and 

= P,{z) , (B12) 

where Pj{z) are the Legendre polynomials. Thus, we 

obtain the familiar three-dimensional expansion 

^ 00 

G^^\ip;x,x') = — ^(2j-f l)P,-(cos6i)Gj(v3;r,r') . 

^ 1=0 

(B13) 

Finally, for d = 4, k = 1, and 

(i i) 2 7 + 2 

(B14) 

VTT (j + 2 ) 3 
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Completing the square in the centrifugal potential in 
Eq. (B5), we make the following approximation for large 
R: 


j{j + d-2) _ {j + nf _ 

4i?2 


i?2 


i?2 


(j + nf 

i?2 


, (B19) 


where k = d/2 — 1, as before. We may then promote 
(i + n)/Rio a. continuous variable k, obtaining 

1 

= — dA:G('u, m), (B20a) 

Jo 

1 

G^^\ip]x,x) = — / dk kG{u,m) , (B20b) 

27r Jo 

G^'^^ipiXjx) = —^ / dkk^G{u,m), (B20c) 
27r Jo 


d-l 


G^'^\lp]x,x) = 


2{Atz)- 2 


dk k'^ ^ G{u, m) , 


(B20d) 


where we have used the general notation employed in 
Sec. Ill, see Eq. (38), with m given by Eq. (40). We 
note that for d > 4, we have dropped terms 0{k/R) and 
higher within the integrand. 


3. Radial function 

For large i?, we neglect the damping term in the ra¬ 
dial equation [Eq. (B5)] and set r = i? in the centrifugal 
potential and discontinuity, giving 


, j'(j +d-2) 


dr 2 


-b 


i?2 


+U"{^) 


Gj{if,r,r') = 


S{r — r') 


(B21) 

Since the solution depends only on the normalized bounce 
u = tanh[ 7 (r — i?)], it is convenient to define 

G{u,u',m) = R‘^~^ Gj{ip;r,r') , (B22) 

cf. Sec. III. Equation (B21) may then be recast in the 
form 




= — 7 ^(I(m — u') , 


G{u, u', m) 

(B23) 


where 


n = 2, = 2 ([l + . (b24) 

Splitting around the discontinuity at u = u', we de¬ 
compose 

G{u,u',m) = ^{u — u')G^ {u,u',m) 

+ i){u'— u)G'^{u,u\m) , (B25) 


where G^{u, u', m) are the solutions to the homogeneous 
equation 




G^{u,u\m) = 0 . 

(B26) 

The latter is the associated Legendre differential equation 
and we obtain the general solutions 

G^(u,m',to) = bl^P™(u)-bB^Q™(u) , (B27) 

where P/J{z) and QJ/iz) are the associated Legendre 
functions of the first and second kind, respectively. 

Matching around the delta function in the inhomoge¬ 
neous equation, we require 

(24> - A<)P^(u') + {B> - B<)Q^{u') = 0 , (B28a) 


+ (b>-b<)To;v)=-:^(54^ 


. (B28b) 


Thus, we find 

= 

B> -B< = 


1 _ QTiu') 2 , 

W[P^{u'),Q^{u')]’ ^ ’ 

_^_ PIM _ (B29b) 

7(i-u' 2) w[g^(u0,d^2™K)] 


where W[P/J{z),Q'!//{z)] is the Wronskian, having the 
explicit form 


w[p;r{z),Q//{z)] = 


{n-m + l) 2 m 
l-u'2 


(B30) 


with the Pochhammer symbol defined in Eq. (B7). We 
also require the boundary condition that G(u,u',m) go 
to zero as u —>■ ± 1, giving 

= — ^cotruTT, P< = 0. (B31) 

We may now solve for the remaining non-zero coeffi¬ 
cients and obtain 


G>{u,u',m) = ^ 


27 sinTOTT 


P2^{u)P^{u') , (B32) 


with G^{u,u',m) = G^{u',u,m). Here, we have used 
the identity 

n{n-m+lh^ P„-”-(^) = ^ cot run P/J (z) - Q//(z) . 
2 sm rmr 2 

(B33) 

Finally, we employ the representation 

m 

Pn{z) = (^^y{n-m+lUPi-^^^\z) (B34) 
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of the associated Legendre function of the first kind in 
terms of the Jacobi polynomials. For n = 2, the polyno¬ 
mial expansion of the latter terminates, and we have 


p^(±m.=Fm) _ 1 (X±m)(2±m) 

— 3(2 ± to)(1 — u) + 3(1 — uY 


(B35) 


After some algebraic simplification, we then arrive at the 
final analytic solution, as presented in Eq. (36) of Sec. III. 


4. Functional determinant 

The normalized heat kernel IC{ip; z, z', k|r), see Sec. Ill, 
is given in terms of the inverse Laplace transform 


K{ip;z,z',k\T) = /:/[G'(M,m)](r) , (B36) 


where 


(B37) 


with 


m = 2 1 + 


-f 
472 


(B38) 


Generalizing to d dimensions, using the continuum 
limit in Eq. (B20d), the correction to the bounce action 
arising from the functional determinant is therefore 






poo 

/ dr 7(1 - 
^0 




n—1 


I -I- u" 


- I^er(”'-4)^erf(n7v^) , (B44) 


where fid = 27r‘^/^/r((i/2) is the solid angle subtended by 
the (d— 1) dimensional hypersphere. The integral over r 
is dominated by r ^ i?, such that (for n = 1,2) 


(-1)"^ drr‘^-7(l-u2) 

We are then left with 

d-l 


1 +u'^ 


-1 


^ r>d—l 

■ 3 ^ ■ 

(B45) 


„(i) _ fld(47r) 2 

r(7i) 

n—1 


Dd—1 /*oo 

- dfc fc'’ 

Jo Jo 


d-2 


' dr 


,-fe r 


erf (nyVr) , 


(B46) 


The inverse Laplace transform may be performed by 
using the shift, scaling and division properties 

C1[F(s + 5)](t) = e-+(r), (B39a) 

CM^M](t) = i/(T/a), (B39b) 

/:"^[s"^E(s)](t) = f dr' /{t') , (B39c) 

Jo 


cf. the form presented in Ref. [86]. 

We may now proceed in one of two ways: (i) perform¬ 
ing the r integration first, we must regularize the k in¬ 
tegral, for instance by introducing an ultraviolet cutoff 
A; or (ii) performing the k integral first, we must instead 
regularize the t integral. The latter is the approach pre¬ 
sented in Ref. [86], which we reproduce in what follows 
for comparison. 


where /(t) = ^[E(s)](t), as well as the elementary 

transformation 

CY^[s-Yir) = ^ , Rez>0. (B40) 
I [z) 

We find 

CY^[m-\m^-n^)-Y{T) = — e[T''(”'-4)-'='l"erf(n7v^) , 

(B4I) 

where 


(i) Performing the r integral first gives 


+ 7*2 -72(n2 - 4) 


(B47) 


Subsequently, performing the k integral for c? = 4, we 
obtain the result in Eq. (79). 


9 

erf(z) = / dt e~* (B42) 

V'^Jo 

is the error function. Hence, we have 

Kiifiz,^,^^) = - ^ 7(1 - 

X - de^'(”'-4)"erf(n7v^) . (B43) 

n—1 ^ ' 


(ii) Instead, performing instead the k integral first, we 
obtain 

1 1 <tzi r°° d±i 

^(47r) 2 / dr r 2 

2 Jo 

2 

X ^ e'^'(”'-4)^erf(n777) , (B48) 

n—1 

which is regularized by introducing a large mass M as 
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follows: 

lim — ^ [ dr 

2 ^ ^ e^oderie)Jo 


X^-^+e^e7^(n=-4). 


erf(n 7 VT) . (B49) 


We may proceed by using the series representation of the 
error function 


erfW = 


( 2 ^+ 1 )!! 


where !! denotes the double factorial. The r integral may 
now be performed, and we obtain 

(B.) 

n=le=0 ^ 

Considering the derivative with respect to e, we have 

d r(e + £+l-d/2)/M2Y 

de fX^) \^) 

r{e + e+l-d/2) /M^Y 

- [^) 

X —'ip{e) + 4’{t + I+l — d/2) , (B52) 

where 'ip{z) is the digamma function. In order to take 
the limit e —)■ 0 safely, we must take note of the poles oc¬ 
curring in r(z) and for non-positive integers. Such 
poles occur in even dimensions for £ = 0, 1, ..., d — S. 

After treating the limit e —>■ 0, we find for d odd (in¬ 
cluding d = 1) 

xEE (LVlV ^(^+1-^/^)- (B53) 

On the other hand, for d even, we hnd 
= -(7i?)''-if2d7r-i 

n=l le=d-2 ^ ' 


where 


= E 


are the harmonic numbers, which we have supplemented 
with iJo = 0 for notational simplicity. 

For d = 4, we then obtain 

2 r °° /9'|2^-|-1 


n=l L^=2 




Lastly, performing the summations, we arrive at the re¬ 
sult 

Defining the counterterms in the proper-time repre¬ 
sentation, see Ref. [ 86 ], and fixing the renormalization 
conditions as in Eq. (64), we find the counterterms 




(B58a) 


(B58b) 


giving 


dB(i) = b( f In ^ - 23 

I 167r2 I \ M2 


Adding these to Eq. (B57), we obtain agreement with 
Eq. (79). 


Appendix C: Renormalization of the AT-field model 

In this final appendix, we highlight the main techni¬ 
cal details of the derivation of the Green’s function and 
corrections to the bounce from the x fields. 

Proceeding as for the isolated ip case, see Sec. Ill, the 
renormalization is fixed using the CW effective poten¬ 
tial [53], evaluated in a homogeneous false vacuum. The 
renormalization conditions are then 


^2Xn/2)2^+i (-l)<i/2-^-i / m2 

( 2 £-fl)!! (d/ 2 -£-I)! V '^472 di2-t-i] , 


d'^UeS 

v=r,x.=o 

a^c/eff 

(9y2 ° 

^ A-i (p—v,Xi—^ 

v=R.x,=o 

d^UeS 

d^UeS 

dxPxj v=v,x.=o 


= 47 " , 


= 67^ -I- ml 
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where the effective potential is 
C/eft = U{ip,x) + ^U{ip,x) 


+ 


with 


^ ^ Ml-kj 

^ ^ dfc +M 2 + ^/c 2 + Ml - 2 k^ , 

(C2) 




^ ' 2 
2 


M2 = - 2 x^ + -^^ + -xi , 


Mi = 


Ml + Ml 


± 


(Ml-Ml 

V 2 


+ aV^x- 


(C3a) 

(C3b) 

1/2 

(C3c) 


and 


ca(v5,x) = + 

+ ^V+4A/22^^ 

^U{(p,x) = + ^Smlxf 

+ ^ ^ 5\^xV^xl + I ^Kxlxl 


(C4a) 


(C4b) 


In Eqs. (C2)-(C4), the summations over i,j = 1,...,N 
have been left implicit for notational convenience. 

Solving the resulting system, we obtain the set of coun¬ 
terterms 

K = 

^ 167r2\72 A2 


-b 


2 I 672 

mi -b 672 mi + 272 


36 O 72 , 67 

In 


2 -L 771-2 


472 


(C5a) 


A 2 A 72 


(w-n) 


/A 2 _ 
V7^ 


30 - In 


r 


/v!!kfin^ 
272 V 4A2 


m: 


-b 1 -b 21N 


A 2 

m 2 -b 272 ^2 
m 2 -b 672 


AX 3A2 

327r2 


In + 5(A^ + 1 ) 


672 -b mi 

jVln ^ - 3IV 


4A2 


f 

X QArl X 


■ 2 r 

\mi + 67 ^ 


2\2 


72 67 - 

InM +41n — 
A 2 


mt 


-b 


13672 , 672 

In 


■ mt 


mi + 272 


472 


4A2 

m 2 — 272 
mi + 672 


(C5b) 


(C5c) 


(C5d) 


A2 /1 72 (m2 - IO72) -b 4327 ^ 

^ “ “3^12^'"A2'^ ^ on2 


(m 2 -b 272 )" 


■ 247 ^ 


(m 2 — 272 )^ — 1127 ^ 672 -b 


(m 2 -b 272 ) ■■ 


In- 


m: 


472 


. (C5e) 


Proceeding as for cp, we find the unrenormalized tad¬ 
pole contribution of the x fields 


E(u) = 


7"A 


16772 
— 3(1 — u^) In 


72 

672 


67 


2 I 2 

m^ 


272 


In 


67 


2 I 2 

m^ 


4A2 


-bl 


■ mt 


4A2 


-6(l-zi2)^(-l)"(n-l-zi2) 

'672 


Jl‘ 2 r ^2 


n—1 

mi \2 / 67 

^ — 1 arccot' 


2 -L m2 \2 
^ - 1 


n 2 j 2 


(C 6 ) 


After adding the counterterms, we obtain 


E«(n) = 


372 A 

167r2 

2 




-2(l-„»)^(-l)"(„-l-,7 

'672 


n 2 j 2 


n—1 

ml \2 f 67 

— 1 arccot 


2 m 2 \2 
^ - 1 


n 2 j 2 


■ (C7) 


We note that the expression in Eq. (C7) agrees with the 
renormalized tadpole contribution from <I> in Eq. (67) for 
m2 = —/i2, as we would expect. Assuming m^ 72 , we 
may expand Eq. (C7) to leading order in giving 

Eq. (83). 

The one-loop correction to the bounce action from 
the determinant over the quadratic fluctuations in the 
X fields is given by 

, , , 9 / ‘ ! 

. 


f\kk^ [ dc^3 

2 io Jo r Jo 


X/:/[S'(u,m)](T) , (C 8 ) 

where S{u,m) is obtained from Eqs. (38) and (75) with 

1 

(C9) 


fc2-bs-bm2\2 

m = vT 1 -b ^ 


672 


Continuing as in Sec. Ill, we find 

■ A2 m2 -b 472 

272 


b9) = 

^ 2 


67 


2 1 2 

' mi 
1 A 2 


mr 


272 


X arccot 


272 3 

^672 




67 ' 


2 ml \2 

^ - I 


n—1 

ml \2 


n2j2 


\ 71,272 


- 1 


(CIO) 
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Adding the counterterm 




^ NR^j' 
2 


3,,3 I 

r 

'x+47^ 

272 


^-20+^ 


21 


mZ 


■2j 


2\2 


In 


67 ^ + m: 
4A2 


672 ^ 

(Cll) 


The result in Eq. (C12) reduces to that found in Eq. (79) 
for = — 27 ^ and N = 1, as we would expect. Instead, 
taking ^ 7 ^, we obtain the expression in Eq. (92). 


obtained in analogy with Eq. (78), we find 


SW = -N 


77373 


ml + 27 ^ 

63-4 , - 63 


272 


mr 


27 ^ 


2\2 


2 ^ 

-V, 

n—1 


67 ^ + m: 

n‘^j2 


2 

X _i 


arccot 


' 67 ^ 


m: 


67 ^ 


n2^2 


- 1 


(C12) 
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